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Nuclear matter properties are calculated in the relativistic mean field theory by using a num- 
ber of different parameter sets. The result shows that the volume energy a± and the symmetry 
energy J are around the acceptable values 16MeV and 30MeV respectively; the incompressibility 
A'o is unacceptably high in the linear model, but assumes reasonable value if nonlinear terms are 
included; the density symmetry L is around lOOMeV for most parameter sets, and the symmetry 
incompressibility K s has positive sign which is opposite to expectations based on the nonrelativistic 
model. In almost all parameter sets there exists a critical point (p c ,S c ), where the minimum and 
the maximum of the equation of state are coincident and the incompressibility equals zero, falling 
into ranges 0.014fm~ 3 < p c < 0.039fm -3 and 0.74 < 5 C < 0.95; for a few parameter sets there is no 
critical point and the pure neutron matter is predicted to be bound. The maximum mass Mns of 
neutron stars is predicted in the range 2.45Mq < Mns < 3.26M0, the corresponding neutron star 
radius Rns is in the range 12.2km< Rns < 15.1km. 

PACS numbers: 21.65.+f, 24.10.Jv, 26.60.+C 



1. Introduction 



Groundstate nuclear matter properties are specified by the nuclear matter equation of state e(pN,5) which is 
simply the energy per nucleon of nuclear matter given as a function of nucleon density pn and relative neutron 
excess S = (p n — p v )j p^. This equation of state is a fundamental quantity in theories of neutron stars and supernova 
explosions, as well as in theories of nucleus-nucleus collisions at energies where nuclear compressibility comes into play 
jjj . The main measured quantities which can provide information about equation of state are the binding energies and 
other data from finite nuclei. As the finite nuclei are in states near the nuclear matter standard state (p^ = po, S = 0), 
which is defined as the equilibrium state of symmetric nuclear matter with minimum energy per nucleon and called 
also the normal state, our actual knowledge of nuclear matter is mainly about nuclear matter at state close to the 
point (po, 0). In this case, the equation of state can be written approximately as Q || 

e(p N ,S) = - ai + i(tf + K s 6*) (EE^El) 2 + f j + £ (£5^1 s*, (l) 
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which is specified by the standard density po, volume energy a%, symmetry energy J, incompressibility Kq, density 
symmetry L and symmetry incompressibility K s . The most interesting quantity for supernova explosion calculation 
is the nuclear incompressibility Kq which dictates the balance between gravity and internal pressure of the stellar 
system, while the most interesting quantities for heavy- ion collision studies are the nuclear incompressibility Kq and 
the symmetry incompressibility K s which influence the side-flow effects and the isotopic distributions of the collisions, 
respectively. 

There is no direct experimental measurement on these quantities. They can be determined only from data fit 
based on some specific nuclear model. Therefore, our actual knowledge about these quantities is essentially model 
dependent. Nowadays the quantities which are known with reasonable precision are ax, J and Kq, being the last two 
still under active investigation. One of the most sophisticated data fit is given by the nonrelativistic Thomas-Fermi 
statistical model of nuclei with Myers- Swiatecki phenomenological nucleon-nucleon interaction ||. It is a fit to 1654 
ground-state masses of nuclei with N, Z > 8, together with a constraint that ensures agreement with measured values 
of the nuclear surface diffuseness, giving the root-mean-square mass deviation equal to 0.655MeV. The data fits based 
on Skyrme nucleon-nucleon interactions give comparable results ||, whereas a model independent but approximate 
data fit also gives ai, Kq, J and L very close to that obtained by the before mentioned data fit g Q. 

As the a-io-p model of the relativistic mean-field theory is used widely to investigate various nuclear phenomena with 
success [Q- [^0|, it is interesting to calculate these nuclear matter quantities within this model by using the available 
parameter sets, to compare with those obtained by the nonrelativistic model. In addition, as these parameters are 
determined by nuclear ground-state properties, it is also interesting to see what the a-uj-p model can predict for the 
nuclear system under extreme conditions of density and asymmetry. In this case, the most interesting quantities 
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are the location e m = e(p m , 6) of the minimum of the equation of state for given asymmetry 5, and the generalized 
incompressibility K m = K{p m , 5) of the nuclear matter at this state [jy] . Another interesting quantity is the maximum 
mass of neutron stars M^s calculated by the equation of state for neutron matter with 6=1. Actually, to predict 
these properties of nuclear matter under extreme conditions is just one of the main goals in developing a relativistic 
mean- field theory |J . 

The purpose of this paper is to make above mentioned calculation in comparing with results obtained by the 
nonrelativistic model. Section II presents the formalism and formulas used in this calculation. Section III addresses a 
numerical analysis on linear a-oj-p model of the relativistic mean-field theory. The standard nuclear matter properties 
calculated from a number of parameter sets are given in Section IV, and the prediction for cold nuclear matter under 
extreme conditions is made in Section V. Section VI gives the summary. Appendix A displays functions F m (x) and 
/ m (cc) which are useful in the analytical expressions as well as in the numerical calculations. The Bjorken-Drell 
convention for four- vector [jl2| and the natural units with K = c = 1 are used. 

2. Formalism 

The a-u)-p model of the relativistic mean-field theory is specified by the following Lagrangian density |J] : 

C = ^{id fi - g^ - g P T-W) - (M - g a( f)]^ 



mlcj?) ~ hlb(g 
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(2) 



where F^ v = d^u) v — d v u>^, B piy = ^b" — c^b^, tjj, <fi, u> and b M are the nucleon, a, to and p meson fields with 
masses M, m a , and m p , respectively, while g a , g u and g p are the respective coupling constants; 6, c and C3 are 
the nonlinear term coefficients, and r are isospin matrices. The nuclear matter equation of state derived from this 
Lagrangian density can be expressed in terms of the nuclear energy density £ as e = £ / pn — M, and 



£ — £k + £ a + £u + £ p , 



(3) 



/if 4t4 



(4) 



£ rr = M 4 



(5) 
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2M 2 (1 + c^ 2 /ml )2 1 C3 4 MX (1 + c^ 2 /m 2 J 4 ' 
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(6) 



C 2 n 2 

tp ~ 2M 2 ' 

where k p and fc„ are the proton and neutron Fermi momenta respectively, 
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and the function F m (x) is defined as (see Appendix A for details): 

F m (x)= [ dxx 2m Vl + x 2 . (10) 
Jo 

The reduced effective nucleon mass £ and thus the field (f> is determined by 

(i -0 + &<^(i -£) 2 + <(i-£) 3 = % 3 ]T Aik/m, (ii) 

7T * — ' 

i—p.n 



and the field ujq by 



CujRn 1 
Mm u 1 + c 3 w 2 /m 2 



Knowing the equation of state, the following formulas for pressure p and generalized incompressibility K [ )lT[ ] can 
be obtained: 

P = -f + P N ~H~ = ^ fc _ ^ M (Pa -£<y+£u>- ^C 3 Wq + £ p , (13) 



' z—p,n 
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£ d PN ~ 3 £[1 + 26C 2 (1 _£) + 3c C 2 (l-£) 2 ] + Q + 3C^ S /M3' 
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In the previous equation, f' m {x) = df m {x)/dx, and the function f m (x) is defined as (see Appendix A for details) 

/ m (z) = / dx^==. (18) 
Jo yl + r 

At the standard state (pojO), the pressure should be zero, 

p(po,0) = 0, (19) 

and 

if = ^(po,0) = 9(p 2 v ^) Q . (20) 
In addition, the following formulas can be derived: 



3 ~ 2 d5 2 



ld 2 e 1 kl C 2 k 3 P 



o 6 + M 2 £q 3tt 2 M 2 : 



(21) 



3 



( d 3 e \ 
\ PN d^df 2 ) 



dp^dS 2 /o 



= J + 2J„ — < 3 
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The subscript in the above formulas stands for the standard state (po,0), and kp is the nucleon Fermi momentum 
of standard nuclear matter which is related to standard density po and nuclear radius constant rg as 
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Formula (|2lj) is well-known in the literature |7J. It is worthwhile to note that, for the linear model with b = c = 
c 3 = 0, the p-meson in the standard state is nonrclcvant to Eqs.(||), (|ll|) and (|l9|). Thus the parameters C a and 
are the same both for model with or without p-meson, since they are determined by standard density po and volume 
energy a%. This point will be discussed more specifically in the next Section. 

3. Determination of C a , and C p in the linear model 

For the linear model, b — c — C3 = 0, the nuclear energy density (|^) in the standard state (poj 0) is simplified as 



M 4 
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and Eq.(jll|) determining the reduced effective nucleon mass £ becomes 



Wo 



'-fi(k F /CoM). 



(28) 



(29) 



In addition, the equilibrium condition (pj| ) and the expression of incompressibility (^0|) are reduced, respectively, to 



2M 4 Q 



/ 2 (W£oM)-^(l-£ ) 



r 2 n 2 
2M 2 



= 0, 



(30) 



Kn = 6 



r 2 k 3 

7T 2 M 2 



+ 3M^f{(k F /^M) 



M 2 Q 



^3-2£ + Qo 



(31) 
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It can be seen from Eqs.(p8[)-(p2[) that the relevant quantities are hp, C 2 , C 2 and M. Note that £0 is determined 
by Eq.(|29|) and po is related to kp by Eq.(B^). Therefore, as the measured nucleon mass can be taken for M, the 
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composite parameters C 2 and C 2 can be determined completely by the value eo = e(po,0), by using Eqs.(|2§|)-(|3C|) 
together with £ — (e + M)p . The procedure is as follows. 

At the stable equilibrium point (po,0), an equation involving eo, kp, £o and C CT can be obtained if Eqs.(|2g|) and 
( ^p| ) are combined to cancel C u . On the other hand, C a can be solved as a function of kp and £o from Eq.([29|). 
Substituting this function of C a into the above-mentioned equation, the following equation involving eo, kp and £o 
can be derived: 

^h{k F /H Q M) + 2/ 2 (fc F /£oA/) = (eo + (33) 

£o can be calculated from this equation, if the location (po> e o) of stable equilibrium point is chosen as input data. 
Having this £o together with po and eo, C a can be calculated from Eq.(p9j), then C w can be determined from Eq.(p8|) or 
(^p|). Finally, the incompressibility Kq can be obtained from Eq. (|3l|) . It can be shown easily that Eq.(33) is identical 



to Eq.(22) of Ref. |9| which is originally given in Ref. |13| . Numerically, po and thus kp can be expressed in terms 
of the nuclear radius constant ro as Eq.(E7[), while eo can be related to the nuclear volume energy coefficient a\ as 
eo = e(po,0) = — a\. The experimentally acceptable values are || 

r fa 1.14fm, a x w 16MeV. (34) 

The numerical calculation shows that, in the ranges 1.05fm< ro < 1.25fm and 15.5MeV< ax < 16.5MeV, the 
effective mass £ ss 0.54 does not depend on the choice of r and ax sensitively. The composite parameters C 2 . and 
C 2 are sensitive to the choice of ro but not of ax- Fig.l shows C 2 and as a function of ro for given ax — 16MeV. 
Fig. 2 gives the nuclear matter incompressibility Kq calculated by Eq. (pl[) as a function of ax for given ro = 1.14fm. 
It is not sensitive to the choice of ro. Furthermore, Fig. 2 shows that Kq is approximately a linear function of ax, in 
agreement with what is obtained in the macroscopic phenomenological approach to the nuclear matter jj^J. 

In case of the Walecka mode l Ea |. C a — 0, other nuclear matter properties J, L and K s can be calculated also 
from these C a and by Eqs.(|21|)-(|2q). The calculated coefficients J, L and K s are almost constant in the range 
15.5MeV< ax < 16.5MeV for given r = 1.14fm, 

J«20MeV, L«70MeV, K s w 88MeV. (35) 

On the other hand, these coefficients depend on the choice of ro weakly, for given Oi. 

In case p-meson is included also in the model, the composite parameter C p can be determined by measured symmetry 
energy J through Eq. (pl|) . The inclusion of p- meson contributes to the symmetry energy with an extra term J p 
(Eq.([23|)) and to the density symmetry L with an extra term 3J P , while keeping the other coefficients oi, Kq and K s 
unchanged. For symmetry incompressibility K s , it can be seen from Eqs.(p6|) and (|lj) that the p-meson contributes 
with a term — 18J P to — 6L and a term 18J P to (l/2)d 2 K/dS 2 \o, and these extra terms cancel each other. 

4. Standard state nuclear matter properties 

There are many parameter sets for the a-u-p model of the relativistic mean field theory in the literature, some of 
them are listed in Table I, where L-W is taken from the pioneering Walecka linear a-cu model |l5|], L-HS from the 
Horowitz-Serot linear a-uj-p model p6[ , LI, L2 and L3 from Lee et al. fL7|| , L-Z, NL-Z and NL-VT from Rufa et al. 
@, NL1 from Reinhard et al. @, NL2 from Fink et al. @, NL3 and NL3-II from Lalazissis et oi_p|, NLB, NLC 
and NLD from Serot §), NL-B1 and NL-B2 from Boussy et al. |§ ||, NL-RA from Rashdan |j7^NL-SH from 



Sharma et al. |25|, TM1 and TM2 from Sugahara and Toki |2q| . Most of them are collected in Reinhard's review [^0| 
In Table I, and are defined, respectively, as 

g 2 = Mbgl, g 3 = cg*. (36) 

It should be noted that some of these parameter sets are given originally in values of Cj instead of g i: i = a, u, p. In 
this case the values of gi given here are calculated from Cj, nrii and M by Eq. ([)[). It should be noted also that our g p 
is only one half of that defined in Ref. || . 

As mainly nuclear matter properties are concerned in the present calculation, the relevant parameters are only C 2 , 
C 2 , C p) while the meson masses m a , and m p are nonrelevant ones, in case of the linear model. However, in 
nonlinear model the meson mass is able in some case to influence the nuclear matter property For example, the u> 
meson mass m w appears in Eq. (|l2]) and thus has effect on nuclear matter property in the nonlinear model via the 
term (w M w M ) 2 . 
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The standard nuclear matter properties related to these parameter sets are shown in Table II, where all quantities 
are given in MeV, except po which is in fm -3 . In the calculation of a%, Kq, J, L and K s , using formulas given in 
Section II and input parameters listed in Table I, Eqs.(O) and ( |l9| ) should be solved simultaneously at first for £o 
and kp at the standard point. The calculation of d 2 K/dS^\o, in Eq.(|26|) of K s , is made numerically, as its analytical 
expression is too complicated to be derived. The simple numerical average among the nonlinear model sets is given 
as the set (NL), and the Myers-Swiatecki's result Q is shown also as the set MS for comparison. 

po and ai give the location of nuclear matter standard state. Most values of po given in the a-co-p model are lower 
than that of Myers-Swiatecki's, the later corresponds to ro = 1.140fm and agrees with that obtained from elastic 
electron scattering and muonic atom spectroscopy measurements [^?J pi] . Most values of a\ given in the a-cu-p model 
are in the reasonable range around 16MeV, except those of LI, L2, L3, LZ and NL2 sets, which seem too large. Since 
a-i is the leading term in the approximate equation of state ([!]), it is the main parameter in any data fit to nuclear 
masses. However, there is a big fluctuation around 16MeV, as can be seen from Table II. 

Kq and J, the next terms to the leading a\ in the approximate equation of state (|l|), are the fine tune in the data fit 
to nuclear masses, as shown in the droplet model of nuclei p9| . It can be seen from Table II that Kq given in the a-u-p 
model is much larger than that of Myers-Swiatecki's, while J is only about 2/3 of Myers-Swiatecki's, for the linear 
a-uj model; J will be increased if the p- meson is added also to the linear <t-lo model, but Kq keeps the same value. 
This is an inherent character of linear a-uo-p model, as has been shown generally in last Section. In this respect, the 
nonlinear terms are needed in order to reduce the nuclear incompressibility Kq, as supported by the calculated results 
listed in Table II. It is worthwhile to note that, even the value of Kq obtained from different nuclear measurements 
and astrophysical observations are spread over a large range from 180 to 800MeV pc| , most expectations based on 
the nonrelativistic model are around 220MeV @. 

Being terms of order higher than Kq and J in the approximate equation of state ([!]), L and K s belong to the 
superfine tune in the data fit to nuclear masses. Even if most values of L given in the a-u-p model seem to be larger 
than the acceptable one, they are still in the reasonable range around lOOMeV. On the other hand, the values of K s 
are all positive whose sign is opposite to most expectations based on the nonrelativistic model 0]. Experimentally, 
K s obtained from the isoscalar giant-monopole resonance energy is between —566 ± 1350 to 34 ± 159MeV M] , 



5. Prediction for cold nuclear matter under extreme conditions 



The stability condition for the state at minimum of equation of state for given asymmetry 8 is 

p(p m ,S) = 0. (37) 

The solution of this equation for given S gives the location of the minimum p m = p m (5). Knowing this location 
Pm(S), the minimum e m = e(p m ,(S) and the generalized incompressibility at this minimum K m (S) = K(p m ,5) can 
be calculated. Furthermore, the critical point of the equation of state (p c , S c ) can be defined as the point where the 
maximum and the minimum are coincident and thus the curvature of e{p, S c ) versus p equals zero. As the generalized 
incompressibility K {p, S) is proportional to this curvature, we have at the critical point 

K m (S c ) = K(p c ,S c ) = 0. (38) 

This equation together with ( |37j ) can be used to obtain the critical point (p c , S c ). 

Table III lists the calculated critical point (p c , S c ), the corresponding effective nucleon mass M* /M, the energy per 
nucleon e m as well as the generalized incompressibility K m at the critical point. In case there is no critical point, the 
corresponding quantities at the minimum point of the pure neutron matter equation of state with S = 1 are listed. p c 
is in fm -3 units, while e m and K m are in MeV units. The values given by the Myers- Swiatecki equation of state lO] 
are also listed in the last row for comparison. It can be seen that there is no critical point for parameter sets LW, LI, 
L2, L3 and NL-B2. In these cases, there is a minimum for the pure neutron matter equation of state and the bound 
neutron matter is predicted. For other parameter sets, the neutron matter is an unbound gas system. The predicted 
critical point (p c , S c ) is in the ranges 0.014fm~ 3 < p c < 0.039fm~ 3 and 0.74 < S c < 0.95, with the corresponding 
effective nucleon mass in the range 0.87 < M* /M < 0.95. 

In addition, the predicted maximum mass M/vs and the corresponding radius Rns of neutron stars, calculated 
by the Oppenheimer-Volkoff equation, using the a-co-p model equation of state of the relativistic mean field theory 
with the above mentioned parameter sets and (5 = 1, are also shown in Table III. The range of the maximum mass is 
2.45M© < Afjvs _ 3.26M Q , and the range of corresponding star radius is 12.2km< Rns < 15.1km. 

Fig. 3 gives some examples of p m {5), where the solid curve from top to bottom in the middle range of 5 corresponds 
to L-W, L-HS, NL-SH, TM1, NLC, and NL1; the dashed curve corresponds to Myers-Swiatecki' result. One source of 
deviation among these curves comes from the difference in the origin of the curves: po = p m (0)- In cases of L-W and 
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LI po is much higher but others are close or lower than that of Myers-Swiatecki's. However, even if all the curves are 
rescaled to the same po, there still exists large deviation among these curves in the middle range of 5. 

Fig. 4 plots some examples of e m versus S, where the solid curve from left to right on the end of the curve corresponds 
to L-HS, NL-SH, TM1, NL1, NLC, and L-W; the dashed curve corresponds to Myers- Swiatecki' result. All curves are 
close each other in the low asymmetry region, but L-W's is significantly lower than others for 5 > 0.2. 

Fig. 5 is the curve K m versus S calculated by same parameter sets as that of Figs. 3 and 4. The solid curve from 
top to bottom in the middle range of S is by L-W, L-HS, NL-SH, TM1, NLC, and NL1, and the dashed curve is by 
Myers-Swiatecki. The diference between these curves is obvious, even if NLC and NLl's are close to each other as 
well as close to Myers-Swiatecki's. 



In summary, the properties of nuclear matter at standard density po with equal neutron and proton densities, 
Pn = Pp, are calculated at first in the relativistic mean field theory with a variety of parameter sets. The result shows 
that the volume energy a\ and symmetry energy J are around the acceptable value 16MeV and 30MeV respectively, 
the incompressibility Kq is reasonable only for nonlinear model while is unacceptably high for linear model, the density 
symmetry L is around lOOMeV for most parameter sets, and the symmetry incompressibility K s has positive value 
whose sign is opposite to most expectations based on the nonrelativistic model. 

Secondly, the calculation shows that for most parameter sets there exists a critical point (jj c , S c ), where the minimum 
and the maximum of the equation of state are coincident and the incompressibility equals zero, and it falls into ranges 
0.014fm -3 < p c < 0.039fm~ 3 and 0.74 < 6 C < 0.95; while for some parameter sets there is no critical point and 
the pure neutron matter is bound. The deviation among results calculated by different parameter sets is discussed. 
The maximum mass of neutron stars is also calculated with results in the range 2.45M Q < Mms < 3.26M Q . It 
is worthwhile to note that a more realistic calculation, by using a nuclear Thomas-Fermi equation of state, gives a 
maximum mass of neutron stars equal to 3.26M© [f^] . The most of observational neutron star masses are between 



As different parameter sets give results which deviate significantly from one another, in order to extract from them 
more reliable predictions for nuclear matter properties, more sophisticated data fit, especially the data fit to larger 
number of nuclear masses and other measured nuclear data is expected for the nonlinear a-w-p model of relativistic 
mean field theory. 



Functions F m (x) and f m (x) defined below are useful in the analytical expressions and numerical calculations of the 
relativistic mean field theory: 



6. Summary 



1.2 - 1.8M, 



APPENDIX A: 




TO > 1, 



(Al) 




TO > 1. 



(A2) 



The following formulas can be obtained: 



F m (x) = f m (x) + f m +i(x), 



(A3) 



f'n+l( x ) 



(A4) 





(A5) 



F' m {x) = {l + x 2 )f m (x), 



(A6) 



f m (x) = xF^ix) - {2m-l)F m _ 1 (x), 



(A7) 
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f m (x) = -xF' m {x) + 2(m + l)F m {x). 
Some examples of F m (x) and f m (x) are: 

- + 2.T 2 )a;v / l + x 2 + ln(y/l + x 2 - x)}, 



h{x) = ^[x^l + x 2 +ln(Vl + x 2 -x)], 



For x <C 1, we have 



h{x) 



F m {x) 



1 - -x 2 )x^/l + x 2 + ln(\/l + x 2 - x) 



x 2m+l x 2m+3 x 2m+5 



+ 



2m +1 2(2m + 3) 8(2™ + 5) 



T 2m+1 T 2m+3 q„2m+5 

/m(x) = 2m +1 ~ 2(2m + 3) + 8(2m + 5) + 



(A8) 

(A9) 
(A10) 
(All) 

(A12) 
(A13) 
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FIG. 1. Composite parameters C% and Ct as a function of ro for given ai = 16MeV, in the linear a-co-p model. 

FIG. 2. The nuclear matter incompressibility K calculated as a function of a\ for given r = 1.14fm, in the linear a-co-p 
model. 

FIG. 3. Some examples of the location p m {8) of the a-co-p model equation of state. The solid curves from top to bottom in 
the middle range of 5 correspond to L-W, L-HS, NL-SH, TM1, NLC, and NL1, respectively. The dashed curve corresponds to 
Myers-Swiatecki's result. 

FIG. 4. Some examples of the deep e m x <5 of the a-co-p model equation of state. The solid curves from left to right on 
the high 8 range correspond to L-HS, NL-SH, TM1, NL1, NLC, and L-W, respectively. The dashed curve corresponds to 
Myers-Swiatecki's result. 

FIG. 5. The curve K m x 8 calculated by the same a-co-p model parameter sets as shown in Figs. 4 and 5. The solid curves 
from top to bottom in the middle range of 8 correspond to L-W, L-HS, NL-SH, TM1, NLC, and NL1, respectively, and the 
dashed curve corresponds to Myers-Swiatecki's result. 



TABLE I. Some parameter sets of the a-co-p model in the relativistic mean-field theory. See text for details. 



Set 


M 










m p 




9<y 








9p 


92 


9'A 


C3 


L-W 


939.0 


550 


.000 


783 


.000 


763. 


9. 


.57269 


11 


,67114 




,00000 


.00000 


.0000 


.0000 


L-HS 


939.0 


520 


.000 


783 


.000 


770. 


10 


.47026 


13. 


.79966 


4 


,03814 


.00000 


.0000 


.0000 


LI 


938.0 


550 


.000 


783 


.000 


763. 


10 


.29990 


12. 


,59990 




,00000 


.00000 


.0000 


.0000 


L2 


938.0 


546 


.940 


780 


.000 


763. 


11 


.39720 


14. 


,24780 




.00000 


.00000 


.0000 


.0000 


L3 


938.0 


492 


.260 


780 


.000 


763. 


10 


.69200 


14 


,87050 




.00000 


.00000 


.0000 


.0000 


L-Z 


938.9 


551 


.310 


780 


.000 


763. 


11 


.19330 


13. 


,82560 


5 


,44415 


.00000 


.0000 


.0000 


NL1 


938.0 


492 


.250 


795 


.359 


763. 


10 


.13770 


13 


,28460 


4, 


,97570 


12.17240 


-36.2646 


.0000 


NL2 


938.0 


504 


.890 


780 


.000 


763. 


9 


.11122 


11. 


,49280 


5, 


,38660 


2.30404 


13.7844 


.0000 


NL3 


939.0 


508 


.194 


782 


.501 


763. 


10 


.21700 


12. 


,86800 


4, 


,47400 


10.43086 


-28.8849 


.0000 


NL3-II 


939.0 


507 


.680 


781 


.869 


763. 


10 


.20200 


12. 


,85400 


4 


,48000 


10.39100 


-28.9390 


.0000 


NLB 


939.0 


510 


.000 


783 


.000 


770. 


9 


.69588 


12. 


,58890 


4, 


,27200 


2.02714 


1.6667 


.0000 


NL-B1 


938.9 


470 


.000 


783 


.000 


770. 


8. 


.75834 


11. 


,80520 


3 


,75195 


7.51446 


-16.8112 


.0000 


NL-B2 


938.9 


485 


.000 


783 


.000 


770. 


9 


,72687 


12 


,89370 


3 


,52938 


9.47080 


-28.1254 


.0000 


NLC 


939.0 


500 


.800 


783 


.000 


770. 


9 


.75244 


12. 


,20370 


4, 


,32984 


12.66960 


-33.3333 


.0000 


NLD 


939.0 


476 


.700 


783 


.000 


770. 


8. 


.26559 


10. 


,86600 


4, 


,49305 


3.79970 


8.3333 


.0000 


NL-RA 


939.0 


515 


.000 


782 


.600 


763. 


9 


.62661 


11 


,90390 


4, 


,52418 


8.06582 


-16.3173 


.0000 


NL-SH 


939.0 


526 


.059 


783 


.000 


763. 


10 


.44400 


12. 


,94500 


4, 


,38300 


6.90990 


-15.8337 


.0000 


NL-VT 


938.9 


483 


.420 


780 


.000 


763. 


9 


.79084 


12. 


,65660 


4, 


.61319 


13.16500 


-38.1282 


.0000 


NL-Z 


938.9 


488 


.670 


780 


.000 


763. 


10 


.05530 


12. 


,90860 


4 


.84944 


13.50720 


-40.2243 


.0000 


TM1 


938.0 


511 


.198 


783 


.000 


770. 


10 


.02890 


12 


,61390 


4, 


.63220 


7.23250 


.6183 


71.3075 


TM2 


938.0 


526 


.443 


783 


.000 


770. 


11 


.46940 


14. 


,63770 


4 


,67830 


4.44400 


4.6076 


84.5318 
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TABLE II. Standard nuclear matter properties , 


given by the 


a-co-p model parameter 


sets listed 


in Table I. oi, K , J, L and 


K s are in 


McV, po in fm 3 . 


See text for details. 










oet 


Po 




-fM) 


T 

J 


T 

L 




T WT 


i 007 

.lyo/ 


1 K 7C 
10. (0 


545.6 


22.11 


7/1 C 

74.5 


7/1 O 

74.8 


T UC 


1 /I o c 

.1450 


1 C 7C 

15. (0 


546.8 


34.98 


115.5 


93.4 


T 1 

Lit 


i 7ftft 
.1 / DO 


i o r.o 
18.02 


625.6 


21.68 


75.6 


81.8 


T O 


1/117 

.141 / 


1 ft 70 

10. (8 


578.5 


19.07 


68.8 


97.4 


T O 


A O A A 

.lo44 


1 O O /I 

18.24 


624.5 


18.86 


69.5 


102.1 


T V 

Li-L 


1 /I O/l 

.14y4 


i 7/ n7 
1 / .U / 


586.3 


48.84 


157.9 


94.2 


\TT 1 


1 ri o 

.1018 


Ad A O 

10.42 


211.1 


43.46 


140.1 


142.6 


MT O 


1 A £ft 
.1400 


i 7/ no. 
1 / .y)6 


onn a 
399.4 


AO Oft 

40.86 


i on 7 

i2y. / 


on i 
2U.1 


MT Q 


1 /I CO 

.1482 


1 ft O /I 

10.24 


271.6 


37.40 


118.5 


100.8 


MT Q TT 


1 /I m 

.i4yi 


1 ft Oft 

10.20 


271.7 


37.70 


119.7 


103.3 


MT T) 


1 /I c c 
.1480 


1 K 77 
10. / / 


a oi n 
421. U 


or ri1 

35.01 


i no o 
1U8.3 


K /I O 

04.8 


MT "P 1 
IN 1j-J3 1 


.1020 


i k 7n 

lo. /y 


280.4 


33.04 


102.5 


76.1 


MT DO 


.102/ 


i £ 7n 

to. /y 


245.6 


33.10 


111.3 


158.8 


mt 


1 /I O K 

.1480 


1 £ 77 
10. / / 


OO A A 

224.4 


ok no 
35.02 


i no n 
108.0 


7ft O 

/0.8 


MT Tl 


.1480 


1 £ 77 
10. / / 


343.2 


35.01 


101.5 


13.5 


MT "D A 

IN I ,-K /\ 


i £7n 
.lO/U 


1 ft O K 

10.20 


320.5 


38.90 


119.1 


62.0 


MT QT-T 


i /iftn 
.140U 


1 ft ok 
10. oO 


o c c o 
355.3 


Oft 1 o 

36.12 


1 1 O ft 

113.6 


7n 7 
/9. / 


NL-VT 


.1530 


16.09 


1 70 O 

1/2.8 


on 70 
3y. / 3 


i oft n 
120.9 


i on n 
130.0 


NL-Z 


.1508 


16.19 


172.8 


41.72 


133.9 


140.0 


TM1 


.1452 


16.26 


281.2 


36.89 


110.8 


33.5 


TM2 


.1323 


16.16 


343.8 


35.98 


113.0 


56.0 


<NL) 


.1500 


16.14 


287.7 


37.53 


117.1 


83.2 


MS 


.1611 


16.24 


234.4 


32.65 


49.9 


-147.1 



TABLE III. Nuclear matter properties at the critical point (p c ,8 c ) or (p m , 1), the maximum neutron star mass Mns and 
the corresponding star radius Rns, calculated by the a-uj-p model parameter sets listed in Table I. p c is in fm~ 3 , e c and K c in 



MeV, Mns 


in solar mass M© 


and Rns in 


km. Myers-Swiatecki's 


values are listed 


in the last 


row for comparison. 


See text for 


details. 
















Set 


S c 


Pc 


M*/M 


e c 


K c 


Mats 


Rns 


L-W 


1.00 


.0987 


.766 


1.93 


77.8 


2.60 


12.2 


L-HS 


.86 


.0392 


.872 


2.75 


0.0 


3.08 


14.6 


LI 


1.00 


.1034 


.718 


-0.63 


142.3 


2.80 


13.0 


L2 


1.00 


.0849 


.712 


-1.04 


138.4 


3.13 


14.4 


L3 


1.00 


.0847 


.688 


-2.46 


174.5 


3.26 


15.0 


L-Z 


.75 


.0388 


.871 


2.51 


0.0 


3.16 


15.1 


NL1 


.91 


.0150 


.951 


1.43 


0.0 


2.96 


14.2 


NL2 


.81 


.0288 


.925 


2.41 


0.0 


2.78 


13.9 


NL3 


.92 


.0182 


.943 


1.62 


0.0 


2.91 


13.9 


NL3-II 


.92 


.0178 


.944 


1.62 


0.0 


2.91 


13.9 


NLB 


.87 


.0327 


.906 


2.43 


0.0 


2.87 


13.8 


NL-B1 


.95 


.0236 


.936 


1.94 


0.0 


2.68 


12.9 


NL-B2 


1.00 


.0212 


.934 


1.75 


1.8 


2.87 


13.5 


NLC 


.95 


.0175 


.949 


1.63 


0.0 


2.77 


13.2 


NLD 


.87 


.0302 


.929 


2.29 


0.0 


2.60 


13.0 


NL-RA 


.87 


.0243 


.934 


1.89 


0.0 


2.75 


13.4 


NL-SH 


.90 


.0235 


.927 


1.90 


0.0 


2.93 


14.1 


NL-VT 


.95 


.0151 


.952 


1.44 


0.0 


2.87 


13.7 


NL-Z 


.94 


.0144 


.953 


1.39 


0.0 


2.92 


13.9 


TM1 


.90 


.0217 


.935 


1.82 


0.0 


2.45 


13.3 


TM2 


.90 


.0217 


.918 


1.83 


0.0 


2.73 


14.4 


MS 


.82 


.0304 




1.10 


0.0 
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500 




100 — 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1.05 1.10 1.15 1.20 1.25 

r (fm) 



0.20 




0.00 — 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
0.0 0.2 0.4 0.6 0.8 1.0 



5 




5 



600 




0.0 0.2 0.4 0.6 0.8 1.0 



5 



